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Noether’s  Theorem 


links  SYMMETRIES  and  conservation  laws 
for  Euler  Lagrange  Systems. 

What  is  a  conservation  law? 

Answer:  a  divergence  expression  which  is 
zero  on  solutions  of  the  system. 

The  heat  equation 

Uf  (  Ux)x  0 

is  its  own  conservation  law.  Integrating, 


where  we  assume  u  is  sufficiently  nice  that  we 
can  interchange  dt  and  f,  and  we  have 
applied  Stokes’  Theorem.  In  words: 


Rate  of  change 


Net  of  comings  and  goings 


of  total  heat  in  Q  across  the  boundary 

no  sources  or  sinks 


The  usual  examples 


Symmetry 
leaves  Ldx  invariant 


Conserved  Quantity 
the  quantity  behind 


the 


D 

d  t 


in  the  Divergence 


t*  —  t  +  c 

translation  in  time 


Energy 


X*  —  t  +  c 

translation  in  space 


Linear  Momenta 


x*  =  Tlx 

rotation  in  space 

a*  —  b) 

b*  —  b ) 

Mb  ~  Ma  =  1 

Particle  relabelling 


Angular  Momenta 


Potential  vorticity 


Variational  Complexes  1-2-3  ! 

are  locally  exact 


SMOG  cf.  P.J.  Olver,  Applications  .  . . 
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Exactness  can  be  used  to  find  conservation 
laws  for  non  Euler-Lagrange  systems  via 
clever  ansatze! 


cf  Hereman,  Sanders,  Sayers  and  Wang,  CRM 
Proceedings ;  Hydon  J.  Phys.  A 


Exactness  is  proved  by  the  use  of  so-called 
homotopy  operators  Hi , 


A2 


Div 


A' 


E 

<r~ 


H  i  H0 


— > 
<r~ 


which  satisfy 


(Divtfi  +H0E)lo  =  oj,  all  w£A3 
Thus  if  E(u>)  =  0,  then  u>  =  Div(iJi(u;)). 


Idea:  solve  ^(clever  ansatz)  =  0  for 
parameters  and  arbitrary  functions.  Then  you 
have  a  conservation  law  using  Hi. 


More  on  d  and  n 


SMOOTH 

d(Ldx)  =  d  Q  (u%  +  d x'j 

—  ( Ux^Ux  +  Uxx^Uxx)dx 

— —  (  'W'CCxd'U  “1“  UxXXX^T)(\x 

H”  (ux^U  —  2.UxxdUx  +  (^ccxdix)^ 
=  £(L)ditda:  +  ^VL 

General  Formula,  explicit,  exact,  symbolic,  for 
r]L  known. 

E  =  7r  o  d,  where  n  projects  out  the 
divergence  term. 

More  than  one  dependent  variable 
6L(x,  u,  v, . .  .)dcc  =  Eu(L)du6xEv (L)6vdx 

^Ux^L 


More  on  d  and  n 


DISCRETE 

d(Lda;)  =  0  +  unun+ij  An 

=  (lindMrt  Un^2^Un  “I-  MndMn_|_2) 
—  (un  +  «n-f2  +  Un_2)d'UnAn 

+(5  —  id)(-  •  •) 


—  -E(Z/n)dunAn  -)-  A(^n) 

General  Formula,  explicit,  exact,  symbolic,  for 
T7£n  known. 


£;  =  7T  o  d,  where  n  projects  out  the  total 
difference  term. 


More  than  one  dependent  variable 

d(f^j iAj-j)  =  Eu(^Lji)(SufiAji  -\-  Ev Aj-j 

+^(r/Ln) 


Variational  Symmetries 


Symmetries  arise  from  Lie  group  actions. 


EXAMPLE:  G  =  (R,  +) 


x 


e  •  x  —  x 


1  —  ex 


_  *  /■  *  \ 

e  •  u  =  u  {x  )  = 


u(x ) 

1  —  ex 


Group  Action  Property 


5  •  (e  •  x)  =  5 


x 


#  \  _  1— ex 


X 


and  similarly  for  u(x). 


\l—exj  1—5 

=  (e  +  5)  •  x 


X 


1— ex 


Prolonged  Group  Action 

du*(x*)  .dx 


* 


Mx 


e  •  Ux  —  tx x *  — 

and 

5  •  (e  •  Ux)  — 


dx  dx  (1  —  ex)2 


5  •  Ux 


^X 


(1  -e(S-x))2  (1  ~(S  +  e)x)2 


Action  on  Integrals 


e  •  Jq  L(x,  n,  ux , . . .)  dx 

def’n  of 

f  ■ 

=  Je.Q  L(e  •  x,e  •  u,e  •  ux ,  •  •  •)  de  •  x 

change  of 
variable 

de  • 

=  Jq  L(e  •  #,  e  •  n,  e  •  ncc,  •  •  •) - dx 

dx 

Use  L2  theory  to  get  that  a  variational 
symmetry  of  a  Lagrangian  is  a  group  action 
such  that 

d  €.  •  X 

L(x,  n,  ux , . . .)  =  L(e  •  x:  e  •  n,  e  •  ux,  •  •  •) - 

dx 


Infinitesimal  Action  on  Integrals 


Since  the  symmetry  invariance  condition 

dc  •  x 

L(x,  u ,  uX: . . .)  =  L{e  •  x,  e  •  n,  e  •  n#,  •  •  •) - 

dx 

is  true  all  e,  then  if  everything  is  sufficiently 
smooth,  applying  — |e=0  to  both  sides,  and 
noting  that  when  e  =  0  we  have  the  identity 
action, 


0 


dL 

dux 


<PX  +  •  •  •  +  L^x 


_  D  |  dL  |  dL  D Q  _i_  i  dL  D^Q  _i_ 

Dec  '  ~du^  '  dux  Dec  '  duxx  D.t2 


where 


Q  =  4>-  ux £,  (p 


de 


e=0e  '  u 


£  = 


de 


e=06 


and  ^  is  the  total  derivative  operator. 


0  =  Div(L£)  +  ^(D  JQa) 


dL 

du°j 


Almost  to  the  punchline 


Let 


VQ  =  E<? 


a 


d 


a 


dua 


Then  the  prolongation  is  defined  by 
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P rvg  =  E  °JQa 


a,J 


du°j 


Note 


n,OL  _ 
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du 
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dx i1  ■  ■  ■  dxp 


Then 


E  D  Q 


dL 
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—  DJua 


dv?j 


prvgjdL 


Recall  that  d  is  one  of  the  two  operators  comprising 
the  Euler  Lagrange  operator,  while  the  left  hand  side  is 
a  divergence  if  Q  is  the  characteristic  of  a  symmetry. 


THE  PUNCHLINE 
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=  a  (L)  +  Div(r/L) 


Q  ■  E(L ) 

=  VQjd(L)  +  Div(prvgjT7i) 

If  Q  is  the  characteristic  of  a  symmetry,  we 
have  that 

*qMl)  =  Div(LO 
and  hence  that 


Q  ■  E(L)  =  Div(something) 


Non-trivial  example 


Semi-geostrophic  equations 


Group 


Invariants 


Equations 


a 


=  4>(a,b) 

6*  =  ^(a,  b) 

h  =  (XaVb  ~  XbVa)-1 
dx  =  h(ybda  —  Dadb) 
dy  =  h(-xbda  +  xadb) 
D tx  =  —j/zDthx  -  jh 
D  ty  =  -filthy  +  jh 


Mb  ~  M 


a 


y 


The  Lagrangian  has  4  arbitary  functions 
which  obey  two  conditions.  The  conserved 
quantity  is  potential  vorticity 


DISCRETE  Almost  Punchline 


This  case  is  easier  than  the  smooth  case. 

•  Since  n  cannot  vary  in  a  smooth  way,  the  “mesh 
variables”  xn  are  treated  as  dependent  variables. 

•  The  group  action  commutes  with  shift: 

e  •  Sj(un )  =  eun+j  =  Sje  •  un 
so  no  prolongation  formulae  are  required. 

For  example, 

Un  Un-\-j 

6  •  Un  —  —  y  €  •  'Un_|_j(  —  — 

1  CXf 'i  1  CXji-^-j 


The  symmetry  condition  is: 
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where  ( )*  =  e  ■  ( ). 


Applying 


de 


e=0 


to  both  sides  of  the  symmetry  condition  yields 


°  =  E 
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Setting 


Qn  = 

then  since 


de 


e=0Xn' 


Qn  = 


de 


* 

e=0  n 


Qxn+k  =  Sk(Qx ),  Qun+k  =  SK(Q%) 


_  c'ki'r-vu 


the  equation  above  can  be  written  as 


0  =  Xq^Lu,  Xq  =  J2  5j(Q“) 
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CK,  J 


dun+J 


DISCRETE  Punchline 


^Ex2^  Ex3  ^  A*  — »■  A2 

UJ  W 

Ln  E(Ln)dun 

=  d(Ln)  +  A(?7j,n) 

Q  ■  E(Ln)  =  XqMLu)  +  A(XQ^Ln) 

Again,  we  get  that  if 

-Xgjd(Ln)  =  0 

then 

Q  •  E(Ln)  =  A(something), 

that  is,  a  total  difference  expression  which  is 
zero  on  solutions  of  the  discrete  Euler 
Lagrange  system. 


Nice  example  T.D.  Lee,  Difference  Equations 
and  Conservation  Laws,  J.  Stat.  Phys.,  46  (1987) 

A  difference  model  for  f(^x2  -  V(x))  dt 


Define 


V(n) 
and  take 


1 

xn — 1 


'Xn 
xn—  1 


V (a:)  da; 


tn  r  n  —  1 


-  V(ri) 


\tn  tn—  l) 


The  group  action  is  t*  =  tn  +  e,  with 
invariant.  The  conserved  quantity  is  thus 
“energy”.  Now,  Qtn  =  1  for  all  n,  and  Q*  =  0. 
The  equations  become 


0  =  Et(Ln) 

0  =  Xgjd  (Ln) 


a 

dtn 

d 


Ln-\~  S 


d 


~dt 


n—  1 


L 


n 


<57 


n 


Ln  ~\~ 


d 

dtr,  _ 


L 


n—  1 


n 


as  Ln  is  a  function  of  ( tn  -  tn- 1)- 


It  is  easy  to  see  in  this  case  that 

0  = (s  - ia)  d;L”) 

is  implied  by  the  two  equations,  to  yield 

+n»)=c 

Note  that  the  energy  in  the  smooth  case  is 

l/2i2  +  V. 

Can  regard  the  EL  eqn  for  the  mesh  variables 
as  an  equation  for  a  variable  mesh. 


INTERLUDE 


If  we  know  the  group  action  for  a  particular 
conservation  law,  we  can  "design  in”  that 
conservation  law  into  a  discretisation  by 
taking  a  Lagrangian  composed  of  invariants. 
These  necesarily  satisfy  vq(I)  =  0  or 
Xq(Iti)  =  0.  The  Fels  and  Olver  formulation 
of  moving  frames  is  particularly  helpful  here: 
a  sample  theorem  is 

Discrete  rotation  invariants  in  z2 

Let  (xm  yn) i  {xmi  ym)  be  two  points  in  the 
plane.  Then 

In,m  =  Xnyn  +  Xmllm •>  Jn,m  —  %nym  Xmlln 

are  rotation  invariants.  Moreover,  any  disrete 
rotation  invariant  is  a  function  of  these. 


Made  up  example 


Suppose 

1  2  !/  x2 

Ln  2  H,n- |-1  xn-\-lVn) 

then 


ZT'a? 

E n 


El  = 


°^n,n+l2/n+l  ^n—l^nUn—l 
—  ^n,n-\-lxn-\-l  Jn—l,nxn—l 


Now, 

Qn  =  (Qn?  Qx)  —  (—Uni  xn)  — 

and  thus 


d 

d<9 


(9=0 


(x 


* 

n? 


Qn  '  Efi  °^n,n+l(  2M2/n-f-  1 

+  ^n-l,n(2M2/n-l  +  %n%n — 1 ) 
—  ^n,n+l^n,n+l  “1“  ^n—  l,n^n—  l,n 
id)(  Jn_i7n/n_i7n) 

gives  the  conserved  quantity. 

Note  that  In,m  — —  Im,n  3nd  Jn,m  — —  Jm,n 


Less  easy  example 


Hereman  et  al . ,  Densities,  Symmetries  and  Recursion 
operators  for  nonlinear  DDEs,  CRM  Proceedings 

The  Toda  lattice  in  polynomial  form  is 

!%  i-'n—  1  vn 

Vn  =  Vn(un  Un_|_i) 

The  scaling  symmetry  is  the  basis  for  the 
ansatz  used  to  obtain  the 
differential-difference  conservation  laws, 
which  are  of  the  form 

~zr~Pn  ~\~  (S'  id)  Jn  =  0 

I— y  L 

for  example 
1  3 

Pn  —  ~un~^~u^vn—  1 —  vn)  ?  Jn  —  un— lun^n—  l~\~v 

sJ 

These  results  use  the  ansatz  plus  homotopy 
operator  method  outlined  earlier. 


5  ro 


Summary  of  the  Pattern 
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Div 
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— » 


A*  — >  A*  — > 


UJ 

Ece  Ea(L)6ua 


=  d(L)  + 


VQ 

XQ 
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Q  •  E(L) 


VQ  \  jdL  + 

XQ  I 


•  the  formula  for  r\L  is  explicit,  exact,  symbolic 


•  the  first  summand  is  a  total  derivative  or 
difference  by  the  symmetry  condition 


OK  let's  try  for  a  Neother’s  Theorem 
for  Finite  Element! 

D.  Arnold,  Beijing  ICM  Plenary  talk 

Given  a  system  of  moments  and  sundry  other 
data,  aka  degrees  of  freedom,  that  yield 
projection  operators  such  that  the  diagram 
commutes: 

0  — >  R  ->•  A0  ->■  A1  — >■  A2  — >■  A3  — >  0 

n0  i  ni  i  n2  i  it  4, 

0  R  -»■  t2  ->■  t3  ->■  0 

all  relative  to  some  triangulation. 

Yields  stability!!  A  Lagrangian  is  composed  of 
wedge  products  of  1-,  2-  and  3-  forms. 

Choose  the  discretisation  of  each  to  be  in  the 
relevant  T%.  Then  commutativity  implies 
conditions  for  Brezzi’s  theorem  to  hold. 


In  one  dimension:  with  en  =  (xn,xn+i),  rio  to 
piecewise  linear,  n x  to  piecewise  constant 
with  moment 


an=  I  u(x)ipn(x)  da; 

Jxn 

Commutativity  of  the  diagram 

u  w^da; 

rio  4-  4-  rii 

u \€n  =  AnX  +  Bn  An  =  /x”+1  v! (x)lf}n(x)  6x 

implies 

An  =  u(x)ll>n(x)]xn+1  ~  (  "+1  u(x)lpn(x)  6x 

Jxn 

Note  that 


1. 


is  required  by  the  projection  property. 


A  finite  element  Lagrangian  is  built  up  of 
wedge  products  of  forms  in  Fq,  F\,  F^,  F3. 
Call  this  resulting  space  JF3.  In  each 
top-dimensional  simplex,  denoted  r,  integrate 
to  get 

i  =  5>r(<4,  •••<*£) 

T 

where  is  the  jth  degree  of  freedom  in  r.  l 

can  also  depend  on  mesh  data. 

Can  now  take  d  which  is  the  variation  with 
respect  to  the  aJT. 


EXAMPLE  In  one  dimension, 


0  -MR  — >■  A0  A  A1  — )■  0 

n0  i  n i  i 

0  — y  M  — y  To  — y  T\  — y  0 
n x  is  to  piecewise  constant  functions  with 
moment  u(n )  =  Jx”+2  u{x)ipn{x)  dx  where 


on  (xn, xnjr2),  while  rig  is  to  piecewise  linear 
functions  with  moments 
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’ X 


n+1 
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u(x)  tix 


n  ^'tz-J-2  '  Xn 

that  is,  an,  an_|_ i  are  used  in  (xn,^n-j-2)i 


u  i — y 


O  '^n 
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n+2  xn 


X  + 


xn-\-l~^~xn-\-2 


xn-\-  2  xn 


an 


xn-\- l~\~xn  | 

xn-\- 2  xn  J 
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Very  simple  example 


J^u^dx  projects  to 

e  r+2  ^n(u)-  dx=E2  ((a2n  ~ a2n+i)2 ) 

n  '’x2n  ^  n  \  %2n-\-2  ®2 n  / 

Then 

at2”  =  4S5PS<Clc,2„-da2„+1) 

_  4  /a2n-a2n+l  _  Q2n-1-Q2n^  d 
\  a'2n+2  ^2n  ^2n-|-l  x2n  J 

+(S  -  id)(something) 

The  discrete  Euler  Lagrange  equation  is  then, 
after  “integration”, 

a2n  ~  a2n+l  _ 

"■  ■  c 


^2n-j-2  ^2n 


Look  now  at  the  "Noether  pattern”  for  the 
Finite  Element  variational  complex 

7TOdo  f  ~ 

— ^  J~2  — ^  J~3  — ^  ^  ^ 

w  w 

Lr  £(Lt)+<K»?l) 

=  dLr 

where  5  is  the  mesh  dependent  coboundary 
operator  (recall  <5(/)(r)  =  /(dr)). 

Step  1:  find  tjl  Step  2:  find  vq 

If  then  vqtj6(Lt)  =  ^(something)  we  will 
have  that 

0  =  Qt  •  E(Lt )  +  <5(something). 


Group  actions  on  moments 


The  clue  is  the  variational  symmetry  group 
action  on  /qL(x,u,  ■  ■  ■)  dx 

Define 


e  ■  fTu(x)i/jr(x)  dec 


=  It  e  ■  u(x)ipT(e  ■  x dec 

Example  Recall  the  projective  action 
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e  ■  x  = 


e  •  u(x)  = 


u(x ) 


1  —  ex  1  —  ex 

Then  the  induced  action  on  the  moments 


OLn  — 
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xn 


X' 
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6x,  pn  = 

Jxn 


xn+ 1  u(x) 


6x 
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is 


e  •  CXn  —  Otrii 


€  '  fin  —  fin  —  € OLn 


In  general  for  this  action, 


e  •  Jxn+1  xmu(x )  da: 


f: 


xn-\-l  xlrl _ u{x) _ da: 

Xn  (1— ex)m  1— ex  (1— exp 


—  rxn+l  xmu(  x)  j 
JXn  (1— ex')m~^~^ 


THINK:  if  you  want  a  coherent  scheme  which 
maps  to  itself  under  this  projective  action, 
and  involves  only  a  finite  amount  of  data, 
then  take  your  moments  to  be 


'xn+ 1  u(x ) 


xn 


X 


m 


6x. 


u(x )  i — y 


m  =  3, 4, ...  N. 


CONCLUSIONS 


The  underlying  algebraic  pattern  of  the 
exact  variational  complexes  provide  a 
framework  for  generalisations  of  Noether's 
Theorem  and  conservation  laws  in 
general. 


Symmetry-adapted  moments  would 
appear  to  be  necessary. 


Next:  formulae  for  r]£JT  where 

d(Lr)  =  E(Lr)  +  5(VLt) 

in  terms  of  the  mesh  dependent 
coboundary  operator. 


